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Abstract 

A  novel  fluorescence  microscopy  technique  for  visualizing  the  transport  of  liquid  water  in  unsaturated  hydrophobic  fibrous  media  has 
been  developed  and  is  applied  to  the  gas  diffusion  layer  of  a  PEM  fuel  cell.  In  the  experiments,  fluorescein  dye  solution  is  pumped  through 
the  fibrous  hydrophobic  gas  diffusion  layer  (GDL)  and  imaged  with  fluorescence  microscopy.  Transient  image  intensity  data  is  correlated 
to  the  liquid  surface  height  and  is  analyzed  and  presented  in  the  form  of  three-dimensional  reconstructions  of  the  time-evolving  gas/liquid 
interface  inside  the  fibrous  structure.  The  high  spatial  resolution  of  the  visualization  can  resolve  the  dynamic  transport  of  liquid  water  through 
distinct  pathways,  which  helps  to  refine  understanding  regarding  liquid  water  transport  mechanism  within  these  porous  layers.  The  physical 
observations  suggest  that  the  water  is  not  transported  via  a  converging  capillary  tree  as  suggested  in  prior  work  and  models.  Rather,  transport 
is  dominated  by  fingering  and  channeling.  Based  on  the  physical  insight  obtained  from  the  experiments,  a  new  water  transport  scheme  is 
proposed  as  the  basis  for  developing  improved  models  for  water  transport  in  hydrophobic  gas  diffusion  layers. 
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1.  Introduction 

The  accumulation  of  liquid  product  water  in  polymer  elec¬ 
trolyte  membrane  (PEM)  fuel  cells  currently  limits  the  per¬ 
formance  of  these  devices.  Improved  water  management  is 
needed  to  achieve  the  increase  in  power  densities  required  for 
commercial  acceptance  of  PEM  fuel  cells  in  transportation 
and  other  demanding  applications.  Liquid  water  accumulates 
in  the  porous  gas  diffusion  layer  (GDL)  of  the  electrodes, 
restricting  the  access  of  reactant  gases  to  catalyst  sites.  In  ad¬ 
dition,  liquid  water  collects  in  the  gas  distribution  channels 
and  manifolding  of  the  fuel  cell  stack.  This  flooding  severely 
hinders  the  upper  power  limit  by  starving  the  supply  of  re¬ 
actant  gases.  Moreover,  liquid  water  flooding  accelerates  a 
number  of  fuel  cell  degradation  mechanisms.  Thus,  signif¬ 
icant  performance  gains,  as  well  as  lifetime  enhancement, 
are  possible  with  a  better  understanding  of  the  liquid  water 
transport  phenomena  present  in  the  fuel  cell. 


*  Corresponding  author.  Tel.:  1  250  721  6034;  fax.:  1  250  721  6323. 
E-mail  address:  ndjilali@uvic.ca  (N.  Djilali). 

0378-7753/$  -  see  front  matter  ©  2005  Elsevier  B.V.  All  rights  reserved, 
doi :  10. 1016/j  .jpowsour. 2005 .03 . 1 99 


In  addition  to  the  growing  number  of  experimental  at¬ 
tempts  to  understand  the  structure  of  the  two-phase  flow  in 
PEM  fuel  cells  [1-12],  there  has  been  many  studies  that 
have  employed  computational  models  [13-16].  Computa¬ 
tional  methods  facilitate  observation  of  liquid  water  trans¬ 
port  in  the  regions  of  PEM  fuel  cells  that  are  inaccessible  to 
in  situ  measurements,  such  as  the  entire  membrane  electrode 
assembly.  However,  the  reliability  of  the  computational  mod¬ 
els  is  still  limited  by  the  lack  of  fundamental  understanding 
of  two-phase  flow  in  fibrous  porous  media.  In  particular,  a  se¬ 
rious  limitation  when  considering  hydrophobic  GDLs  is  the 
relationship  used  to  model  capillary  flow.  Typically,  the  cap¬ 
illary  characteristics  are  modeled  with  a  Leverett  J-function 
in  a  form  obtained  experimentally  by  Udell  [17]  for  soil¬ 
like  porous  media.  In  a  recent  review,  Litster  and  Djilali  [18] 
showed  that  there  is  a  significant  difference  between  the  cap¬ 
illary  relationship  of  Udell  and  that  found  for  a  Teflon  coated 
porous  media  similar  to  a  GDL.  Thus,  experimental  charac¬ 
terization  of  capillary  flow  in  GDLs  is  required  before  the 
accuracy  of  two-phase  PEM  fuel  cell  models  can  be  veri¬ 
fied.  A  brief  review  of  GDL  characteristics  and  of  previous 
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visualization  studies  is  presented  below  prior  to  introducing 
the  fluorescence-based  technique  developed  for  this  study. 

1.1.  Gas  diffusion  layers 

PEMFC  electrodes  feature  two  regions  of  porous  media: 
the  gas  diffusion  layer  and  the  comparatively  more  dense  cat¬ 
alyst  layer.  The  gas  diffusion  layer  is  much  thicker  and  more 
porous  than  the  catalyst  layer.  The  catalyst  layer  features  sig¬ 
nificantly  less  void  space  because  of  the  high  content  of  re¬ 
quired  proton  conducting  ionomer  (typically  Nation)  [19]. 
Further  variation  exists  in  the  properties  of  the  gas  diffusion 
layers  employed  in  PEM  fuel  cells.  Two  structures,  both  con¬ 
sisting  of  carbon  fibers,  dominate  the  forms  of  gas  diffusion 
layers  currently  available:  carbon  cloth  and  carbon  paper. 
Carbon  cloth  is  constructed  of  woven  tows  consisting  of  in¬ 
dividual  carbon  fibers.  Alternatively,  carbon  paper  is  formed 
from  randomly  laced  carbon  fibers. 

Carbon  paper  electrodes  feature  a  much  more  porous 
structure  than  their  cloth  counterparts.  The  mean  pore  di¬ 
ameter  and  porosity  of  carbon  paper  is  approximately  35  jam 
and  85%,  respectively.  Carbon  cloth  values  are  much  lower 
and  difficult  to  quantify  [20,21].  Both  Mathias  et  al.  [20]  and 
Ihonen  et  al.  [21]  present  detailed  characterizations  of  gas  dif¬ 
fusion  layer  materials.  Another  issue  with  carbon  cloth  is  that 
it  is  only  available  in  thicknesses  between  350  and  1000  p,m, 
whereas  carbon  paper  is  available  in  thicknesses  as  low  as 
90  (Jim.  Additional  thickness  results  in  greater  power  losses 
due  to  Joule  heating  and  increased  resistance  to  mass  trans¬ 
fer.  In  addition,  the  two  gas  diffusion  layer  structures  vary  by 
spatial  uniformity  and  degree  of  anisotropy.  Carbon  cloth,  be¬ 
cause  of  its  woven  structure,  is  spatially  heterogeneous  on  a 
macroscopic  scale,  whereas  carbon  paper  is  roughly  spatially 
homogeneous  due  to  its  random  lacing.  Moreover,  the  woven 
nature  of  carbon  cloth  results  in  three  degrees  of  macroscopic 
anisotropy,  whereas  carbon  paper  results  in  only  two  degrees. 
All  three  forms  of  porous  media  in  PEM  electrodes  are  sum¬ 
marized  in  Table  1 . 

Generally,  gas  diffusion  layers  are  treated  with  a  PTFE  so¬ 
lution  (Teflon)  to  increase  the  hydrophobicity  of  the  medium. 
This  aids  water  management  in  the  electrode.  The  hydropho¬ 
bicity  causes  water  droplets  to  agglomerate  and  keeps  liquid 
water  from  reentering  the  GDL  after  being  expelled.  Also, 
the  hydrophobicity  reduces  condensation.  Condensation  oc¬ 
curs  in  PEM  fuel  cells  because  the  fuel  and  oxidant  gases  are 
typically  saturated  with  water  vapor.  It  has  been  postulated 
that  the  condensation  takes  place  in  hydrophilic  cracks  in  the 
carbon  fibers  [2] . 


1.2.  Visualization  of  two -phase  flow  in  PEMFC 

Direct  visualization  of  the  transport  of  liquid  water 
through  the  microscale  pores  of  unsaturated  gas  diffusion 
layers  has  yet  to  be  reported.  This  is  largely  due  to  the  im¬ 
mense  difficulty  associated  with  making  in  situ  measure¬ 
ments  during  the  operation  of  PEM  fuel  cells.  This  difficulty 
is  compounded  by  the  requirements  for  microscale  flow  visu¬ 
alization,  which  include  optical  accessibility.  However,  there 
has  been  some  progress  using  three  methods: 

1 .  Ex  situ  visualization  of  water  droplets  placed  or  condensed 
on  an  isolated  GDL  without  net  transport  of  the  liquid  [1- 
4]. 

2.  Use  of  optically  accessible  flow  fields  to  visualize  the 
emergence  of  macrodroplets  of  liquid  water  from  the  sur¬ 
face  of  the  GDL  as  well  as  their  transport  through  the  flow 
channels  [1,5-7]. 

3.  Implementation  of  radiography  techniques  to  macroscop- 
ically  visualize  water  distribution  within  an  operational 
fuel  cell  [8-12]. 

The  first  approach  furthers  understanding  of  the  interac¬ 
tion  between  liquid  water  and  GDLs.  Specifically,  the  level  of 
hydrophobicity  can  be  determined.  However,  external  visu¬ 
alizations  of  static  droplets  do  not  resolve  the  bulk  transport 
of  liquid  water  within  the  porous  structure.  The  last  two  tech¬ 
niques  target  the  transport  of  liquid  water  within  the  fuel  cell 
flow  field.  Nevertheless,  these  techniques  do  offer  the  oppor¬ 
tunity  of  identifying  locations  where  water  is  expelled  from 
the  GDL.  The  critical  nature  of  the  problem  and  the  atten¬ 
dant  challenges  are  highlighted  by  the  fact  that  all  the  studies 
referenced  above  were  published  in  the  last  2  years. 

1.2.1.  Ex  situ  visualization 

To  visualize  the  contact  angle  between  liquid  water  and  gas 
diffusion  layers,  Tuber  et  al.  [1]  and  Pasaogullari  and  Wang 
[3]  captured  images  of  liquid  droplets  resting  on  the  surfaces 
of  GDLs  with  varying  hydrophobicity.  Lim  and  Wang  [4] 
conducted  a  study  of  hydrophobic  polymer  content  in  GDLs, 
quantifying  the  hydrophobicity  with  an  adapted  capillary  rise 
method.  They  also  introduced  an  additional  method  of  visual¬ 
izing  the  contact  angle.  By  utilizing  a  microscope  and  a  CCD 
camera,  Pasaogullari  and  Wang  visualized  the  rise  height  on 
the  exterior  surface  of  a  GDL  dipped  into  a  bath  of  heated 
water. 

Nam  and  Kaviany  achieved,  with  the  aid  of  environmental 
scanning  electron  microscopy  (ESEM),  the  closest  realiza¬ 
tion  of  microscale  liquid  water  transport  visualization  in  an 


Table  1 


Summary  of  porous  media  in  PEM  electrodes  [20,21] 


Porous  media 

Spatial  uniformity 

Dimension  of  anisotropy 

Porosity  (%) 

Mean  pore  diameter  (/xm) 

Thickness  (/xm) 

Catalyst  layer 

Homogeneous 

Isotropic 

5-55 

~  1 

~  15 

Carbon  cloth 

Heterogeneous 

3D 

30-40 

~10 

350-1000 

Carbon  paper 

Homogeneous 

2D 

70-90 

~ 3 0-40 

90-350 
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unsaturated  GDL.  Their  imaging  presented  the  distribution 
of  liquid  droplets  that  had  condensed  on  the  fibers  of  a  hy¬ 
drophobic  GDL.  ESEM  is  a  type  of  SEM  that  places  the  sam¬ 
ple  in  a  low-pressure  condensing  environment  such  as  water 
vapor.  Condensation  is  induced  by  incrementally  increasing 
pressure.  Though  this  visualization  offers  insight  into  the  dis¬ 
tribution  of  condensation  sites  within  a  GDL,  it  does  not  pro¬ 
vide  visualization  of  the  bulk  transport  of  liquid  water  through 
the  GDL.  This  study  provided  nonetheless  some  important 
insight  into  the  local  motion  of  the  water  microdroplets  and 
their  coalescence  inside  hydrophobic  fibrous  media. 

7.2.2.  Optically  accessible  flow  fields 

Optically  accessible  flow  fields  allow  fuel  cell  researchers 
to  visualize  the  transport  of  liquid  water  through  the  gas  de¬ 
livery  channels.  These  techniques  rely  on  the  use  of  CCD 
cameras  to  record  phase  distribution.  Water  transport  visu¬ 
alization  in  flow  fields  has  been  reported  by  the  research 
groups  at  Penn  State  [5,7],  and  at  the  Fraunhofer  Institute 
[1,6].  This  form  of  visualization  can  aid  in  the  understanding 
of  liquid  flows  in  unsaturated  GDLs  by  resolving  the  sub¬ 
tleties  of  droplet  formation  on  the  surface  of  the  GDL  and 
the  distribution  of  emergence  sites.  However,  these  visual¬ 
izations  provide  little  insight  into  the  behavior  or  distribution 
of  liquid  water  within  the  GDL. 

1.2.3.  Radiography 

By  constructing  a  fuel  cell  out  of  non-magnetic  materials, 
Tsushima  et  al.  [8]  used  MRI  to  visualize  the  transport  of 
liquid  water  across  a  polymer  electrolyte  membrane  during 
operation.  The  nature  of  the  MRI  technique  restricted  tests  to 
a  wider  than  average  membrane.  However,  this  did  provide 
insight  into  the  otherwise  invisible  transport  of  water  at  the 
gas  diffusion  electrode/membrane  interface.  Feindel  et  al.  [9] 
recently  overcame  the  limitation  of  visualizing  liquid  water 
in  the  flow  field  with  MRI  by  fabricating  flow  fields  from 
machined  Delrin  with  gold  contacts.  But  because  graphite 
GDLs  were  still  required,  the  water  content  in  the  GDL  could 
not  be  adequately  visualized. 

Neutron  imaging  has  shown  potential  for  visual  inspec¬ 
tion  of  transport  of  liquid  water  in  an  operating  fuel  cell.  The 
relatively  higher  sensitivity  of  neutron  transmission  to  the  hy¬ 
drogen  in  liquid  water  than  to  metals  makes  it  ideal  for  mon¬ 
itoring  two-phase  flows.  Material  limitations  are  few  relative 
to  other  tomography  and  radiography  techniques.  The  use 
of  neutron  imaging  for  visualizing  liquid  water  transport  in 
fuel  cells  was  first  reported  by  Satija  et  al.  [10].  They  utilized 
the  National  Institute  of  Standards  and  Technology’s  (NIST) 
neutron  beam  line  and  were  able  to  generate  videos  of  water 
transport  by  capturing  1000  images  at  2  s  intervals.  However, 
because  of  the  nature  of  radiography,  only  planar  integrals 
of  water  content  were  available.  This  feature  makes  visual¬ 
ization  of  through-plane  flow  in  individual  pores  in  a  GDL 
very  difficult.  Satija  et  al.  used  tomography  to  reconstruct  the 
three-dimensional  structure  of  a  dry  cell  not  in  operation.  A 
significant  disadvantage  of  tomography  for  flow  visualiza¬ 


tion  of  water  transport  in  GDLs  is  the  temporal  resolution, 
which  is  on  the  order  of  seconds  for  good  quality  images. 

Recently,  Kramer  et  al.  [11]  published  their  findings  em¬ 
ploying  neutron  radiography  to  perform  PEM  fuel  cell  di¬ 
agnostics.  Their  paper  presents  a  number  of  advancements 
in  data  treatment  and  water  quantification.  Pekula  et  al.  [12] 
have  also  recently  presented  results  of  an  experimental  analy¬ 
sis  of  water  transport  in  an  operating  PEM  fuel  cell.  Notably, 
they  were  able  to  acquire  images  at  a  rate  of  30  frames  s_1, 
which  is  a  significant  improvement  over  the  values  reported 
by  Satija  et  al.  and  Kramer  et  al.  The  effect  of  this  higher 
frame  rate  on  image  quality  was  not  discussed.  Neverthe¬ 
less,  Pekula  et  al.  were  able  to  identify  water  entrapment  in 
recirculation  zones  located  in  the  bends  of  the  flow  field. 

There  are  three  distinct  reasons  why  most  visualization 
techniques  so  far  have  been  unable  to  achieve  visualization 
of  two-phase  flow  within  GDLs. 

1.  GDLs  are  optically  opaque. 

2.  GDLs  are  formed  of  graphite. 

3.  GDLs  have  characteristic  geometries  with  length  scales 

under  50  pan. 

The  opacity  of  GDLs  eliminates  methods  necessitating 
complete  optical  access,  such  as  conventional  microparticle 
image  velocimetry  (PIV)  ,  photobleaching,  caged  fluores¬ 
cence,  and  other  velocimetry  techniques  that  track  particles 
or  discrete  fluid  volumes  [22].  The  second  issue  is  that  the 
presence  of  graphite  limits  the  use  of  techniques  that  employ 
X-ray  transmission  [10]  such  as  X-ray  micro-PIV  [23]).  The 
third  point  rules  out  the  use  of  transmission  radiography  tech¬ 
niques  due  to  typical  resolutions  greater  than  the  mean  GDL 
pore  diameters  of  35  pan.  In  addition,  radiography  presents 
difficulties  in  resolving  flow  within  a  single  pore  because  of 
the  integral  image.  It  is  likely  that  with  modest  improvements 
in  resolution  the  planar  neutron  imaging  technique  of  Satija 
et  al.  [10]  would  be  very  valuable  for  GDL  analysis.  The 
complexity  and  cost  of  neutron  imaging,  however,  present 
additional  challenges. 

In  this  paper  we  present  a  novel  method  for  visualizing  the 
transport  of  liquid  water  in  unsaturated  gas  diffusion  layers. 
The  visualization  is  sought  to  improve  computational  mod¬ 
els  and  provide  a  new  means  to  evaluate  GDLs.  To  achieve 
these  goals,  fluorescence  microscopy  is  employed  to  visual¬ 
ize  fluorescein  dye  (water-based)  pumped  through  the  GDL 
structure.  This  method  overcomes  the  majority  of  obstacles 
faced  by  other  methods,  but  has  some  limitations.  Through 
image  processing,  three-dimensional  renderings  of  the  tem¬ 
poral  evolution  of  the  gas/liquid  interface  are  constructed. 
The  experimental  observations  are  used  to  postulate  a  new 
water  transport  mechanism  in  hydrophobic  fibrous  GDLs. 

2.  Experimental 

Ex  situ  visualization  of  liquid  water  transport  in  gas  dif¬ 
fusion  layers  is  accomplished  with  fluorescence  microscopy. 
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Fluorescence  imaging  is  a  form  of  scalar-based  flow  visu¬ 
alization  in  which  the  properties  of  a  flow  field  are  eluci¬ 
dated  with  a  conserved  scalar.  Fluorescence  imaging  includes 
a  multitude  of  techniques  employing  certain  molecules  that 
emit  a  photon  immediately  after  being  excited  by  electromag¬ 
netic  radiation.  Poly  aromatic  hydrocarbons  are  commonly 
employed  for  fluorescence  imaging  [22].  The  fluorescence 
process  can  be  idealized  in  three  steps:  absorbance,  dissipa¬ 
tion,  and  emission.  First,  the  molecule  absorbs  a  photon  from 
an  excitation  beam  (blue  light,  A  =  488  nm).  Second,  the  ex¬ 
cited  state  achieved  from  absorbing  the  photon  dissipates  to 
the  surrounding  environment.  Third,  a  photon  of  lower  en¬ 
ergy  is  emitted  (green  light,  A  >  510nm).  Quantitative  fluo¬ 
rescence  imaging  is  a  well-developed  field.  For  example,  the 
local  intensity  of  the  emitted  green  light  can  be  correlated  to 
the  concentration  of  fluorescent  dye  to  quantify  mixing  [24]. 

A  description  follows  of  the  gas  diffusion  layer,  the  exper¬ 
imental  apparatus  and  the  image  acquisition  and  processing 
techniques  used  in  this  study. 

2.7.  Gas  diffusion  layer 

The  present  analysis  employed  a  Toray  TGP-H-060  non- 
woven  fibrous  gas  diffusion  layer.  The  GDL  is  coated  with 
PTFE  for  the  desired  hydrophobic  surface.  The  PTFE  con¬ 
tent  of  the  GDL  studied  was  10  wt.%.  The  TGP-H-060  GDL 
features  a  thickness  of  190  [xm.  The  mean  pore  diameter  of 
this  layer  has  been  determined  to  be  approximately  40  ixm 
from  capillary  flow  porometry  and  mercury  porosimetry  [20]. 
Moreover,  a  mean  fiber  diameter  of  8  [xm  was  been  deter- 


Fig.  1 .  Optical  microscope  image  of  a  Toray  TGP-H-060  gas  diffusion  layer 
with  a  PTFE  content  of  10  wt.%. 

mined  from  microscope  images.  An  image  of  this  layer,  cap¬ 
tured  with  an  optical  microscope,  is  presented  in  Fig.  1. 

2.2.  Apparatus 

The  experimental  apparatus  consisted  of  an  optical  micro¬ 
scope  arranged  for  fluorescence  microscopy,  a  cooled  CCD 
camera,  a  numerically-controlled  syringe  pump,  and  an  as¬ 
sembly  for  pumping  fluorescein  dye  through  the  GDL.  A 
schematic  of  the  experimental  apparatus  is  shown  in  Fig.  2. 
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Fig.  2.  Schematic  of  the  experimental  apparatus. 
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Fig.  3.  GDL  holding  assembly. 


The  GDL  is  held  in  place  above  the  fluid  reservoir  by  the 
compression  of  a  screwed-on  top  plate.  The  top  plate  also 
effectively  seals  the  edge  of  the  GDL.  In  order  to  concentrate 
water  transport  at  the  center  of  the  GDL,  the  bottom  surface 
was  masked  except  for  a  ^  2  mm  x  2  mm  window  in  the  cen¬ 
ter  of  the  apparatus.  A  photograph  of  the  GDL  assembly  is 
shown  in  Fig.  3. 

Images  were  captured  with  an  upright  Leica  Microsys¬ 
tems  DM  LM  optical  microscope  featuring  a  cooled  CCD 
camera.  Magnification  was  achieved  using  an  objective  with 
5x  magnification  and  a  numerical  aperture  (NA)  of  0.15. 
These  characteristics  result  in  a  depth  of  field  (Jfieid)  equal  to 
23  pin. 


X  510nm 

Afield  2  r\  i  c2  L  IP  • 

NA  0.15z 


The  depth  of  field  is  an  important  parameter  in  the  present 
study  because  it  corresponds  to  the  though-plane  height  of  the 
gas  diffusion  layer  that  can  be  in  focus  simultaneously.  With 
respect  to  the  fluorescent  dye,  the  depth  of  field  is  a  measure 
of  the  vertical  extent  of  the  imaging  volume  in  which  the  flu¬ 
orescence  signal  is  integrated.  Since  the  numerical  aperture 
of  microscope  objectives  generally  increases  with  magnifica¬ 
tion,  greater  lateral  resolution  arrives  at  the  cost  of  a  reduced 
depth  of  field.  For  this  reason,  a  low  magnification  of  5  x  was 
chosen  to  maximize  depth  of  field,  while  maintaining  ade¬ 
quate  lateral  resolution.  The  trade-off  between  the  range  of 
the  fluid  height  measurement  and  the  planar  resolution  of  the 
area  in  focus  is  a  limitation  of  this  method. 


2.3.  Image  acquisition  and  processing 


To  visualize  the  liquid  transport,  a  solution  of  water  and 
fluorescein  is  pumped  through  the  plane  of  the  GDL  and  im¬ 


aged  with  fluorescence  microscopy.  The  fluorescent  dye  em¬ 
ployed  here  was  fluorescein,  332  MW  (Molecular  Probes). 
The  dye  was  dissolved  to  a  concentration  of  0.5  mM  in  pure 
water.  Due  to  the  imaging  intensity  of  the  fluorophores  this 
relatively  low  concentration  provided  ample  signal  and  con¬ 
tributed  negligibly  to  the  physical  properties  of  the  solution. 
With  respect  to  surface  properties,  by  imaging  droplets  on 
the  GDL  surface  it  was  determined  that  both  pure  water  and 
water  with  the  fluorescein  marker  exhibited  similar  wetting 
behavior. 

The  path  of  the  liquid  fluorescein  dye  evolves  as  the  pres¬ 
sure  in  the  fluid  reservoir  increases.  The  rate  of  pressure  in¬ 
crease  is  minimized  and  controlled  by  the  pumping  rate  via 
a  syringe  pump.  Because  of  the  slow  pressure  increase,  ob¬ 
servations  of  the  step-by-step  evolution  of  the  fluid  path  are 
made.  Images  of  the  fluorescein  dye  are  captured  at  an  ad¬ 
equate  time  interval  (0.25  s)  to  reveal  the  various  paths  of 
the  water  inside  the  GDL.  Subsequent  image  analysis  cor¬ 
relates  the  light  intensity  to  the  relative  height  of  the  liquid 
water  surface.  Additional  image  processing  provides  three- 
dimensional  renditions  of  the  liquid  water  surface  and  the 
surface  velocities. 

In  many  studies,  such  as  that  by  Oddy  et  al.  [24],  the  flu¬ 
orescence  intensity  field  imaged  is  correlated  to  the  concen¬ 
tration  field  of  the  fluorescent  marker.  In  the  present  study, 
the  presence  of  fluorescence  signal  differentiates  the  liquid 
phase  from  the  gas  phase.  Since  the  concentration  of  fluores¬ 
cent  dye  in  the  liquid  is  uniform,  the  local  relative  intensity 
of  the  fluorescence  signal  imaged  is  assumed  to  vary  pro¬ 
portionally  to  the  depth  of  excited  fluorescein.  Herein,  the 
two-dimensional  intensity  array  recorded  by  the  CCD  cam¬ 
era  is  linearly  correlated  to  the  relative  height  of  the  gas/liquid 
interface. 

The  correlation  between  the  range  of  the  height  and  the 
range  of  the  image  intensity  is  derived  from  the  depth  of 
field  and  experimental  observations.  The  observable  range  of 
the  surface  height  presented  herein  is  30  |xm ,  slightly  greater 
than  the  nominal  23  jxm  depth  of  field.  This  estimated  range 
is  supported  by  observations  of  the  fluorescein  dye  within 
the  height  of  three  to  four  8  |xm  diameter  fibers.  The  linear 
correlation  between  surface  height  ( q )  and  the  8-bit  intensity 
data  (/)  is  expressed  as: 


30  x  10  6  m 
- 1 

255 


Prior  to  correlating  the  image  intensities  to  heights,  im¬ 
ages  were  processed  with  a  smoothing  filter  to  reduce  high 
frequency  noise  in  the  image.  The  signal-to-noise  ratio  (SNR) 
of  this  method  is  relatively  high  when  little  fluorescein  is 
present.  However,  as  more  fluorescein  appears,  the  reflectiv¬ 
ity  of  the  carbon  fibers  becomes  significant  and  decreases 
the  SNR  to  some  extent.  The  contribution  of  reflected  light 
to  the  signal  was  effectively  removed  by  thresholding.  The 
digital  resolution  of  the  acquired  images  is  325  x  258  pixels 
resulting  in  a  spatial  resolution  of  5.35  jxm pixel-1. 
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3.  Results  and  discussion 

Time- sequenced  images  of  the  fluorescein  dye  travel¬ 
ing  through  the  GDL  are  presented  in  Figs.  4  and  5.  The 
images  are  those  captured  directly  using  fluorescence  mi¬ 
croscopy  with  the  image  contrast  inverted  for  clarity.  A  three- 
dimensional  rendering  of  the  image  intensity  data  is  displayed 
next  to  each  captured  image.  The  three-dimensional  render¬ 
ings  feature  a  magnified  vertical  scale  to  improve  the  visibility 


of  surface  features.  The  images  and  renderings  are  presented 
at  intervals  of  14.5  s. 

In  the  images  at  the  start  of  the  sequence  there  is  a  very 
low  disperse  signal  emanating  from  within  the  GDL.  Never¬ 
theless,  the  largest  signal  in  the  first  two  slides  is  located  at 
the  center,  indicating  the  location  of  the  liquid  surface.  On  the 
third  slide,  however,  there  is  a  sudden  appearance  of  floures- 
cein  dye  closer  to  the  top  edge  of  the  images.  Subsequent 
images  and  three-dimensional  renderings  indicate  that  the 
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Fig.  4.  Time  sequenced  digital  images  of  dye  transport  in  the  GDL  and  corresponding  three-dimensional  renderings.  Contrast  of  digital  images  is  inverted. 
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Fig.  5.  Continuation  of  the  time  sequenced  digital  images  of  dye  transport  in  the  GDL  and  corresponding  three-dimensional  representations  from 
Fig.  4. 


surface  of  the  fluid  reached  a  particularly  small  fiber  cross- 
section  and  its  further  progression  is  restricted  by  the  strong 
capillary  pressure. 

After  the  initial  advance  of  water,  a  number  of  additional 
paths  develop  more  slowly  in  the  upward  direction.  It  be¬ 
comes  evident  in  the  images  and  renderings  later  in  the  se¬ 
quence  that  there  is  one  dominant  fluid  path  located  below  the 
path  that  first  emerged.  At  times  between  159.5  and  202.5  s, 
the  liquid  rises  through  this  dominant  path,  while  the  fluid  in 


an  adjacent  path  falls.  This  is  attributed  to  a  local  decrease  in 
pressure  as  the  fluid  in  the  dominant  path  rises  quickly  having 
broken  through  a  local  constriction. 

In  the  final  images  and  three-dimensional  renderings  it 
becomes  evident  that  the  liquid  is  on  the  verge  of  breaking 
through  the  GDL.  In  the  last  image,  it  appears  the  two  top  fluid 
paths  have  merged  into  one  large  path.  This  is  because  the  two 
paths  have  merged  at  the  surface  of  the  GDL.  Immediately 
after  the  last  image  in  Fig.  5,  a  spherical  droplet  forms  on  the 
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Fig.  6.  Raw  fluorescence  microscopy  image  of  a  droplet  immediately  after 
breakthrough. 

surface  the  GDL.  A  raw  fluorescence  microscopy  image  of 
this  droplet  is  shown  in  Fig.  6.  This  figure  clearly  shows  the 
reflectivity  of  the  fibers. 

3.1.  Analysis  of  distinct  flow  paths 

To  further  identify  the  dynamic  interaction  between  adja¬ 
cent  flow  paths,  five  of  them  have  been  selected  for  further 
analysis.  Fig.  7,  which  is  a  contour  plot  of  filter  image  inten¬ 
sity  at  a  time  of  175  s,  presents  the  five  flow  paths  that  have 
been  chosen.  The  paths  discussed  previously  are  paths  1-3. 
Path  1  is  the  route  of  the  first  liquid  advance.  Path  2  is  the 
dominant  pathway,  where  the  liquid  breaks  out  of  the  GDL. 
Path  3  is  the  pathway,  where  the  liquid  surface  is  shown  to 
fall,  while  the  surface  in  path  2  quickly  rises  (see  Fig.  5). 

Fig.  8  presents  the  surface  height  of  each  distinct  path 
versus  time.  This  plot  features  a  much  higher  temporal  reso¬ 
lution  than  that  of  the  snapshots  in  Figs.  4  and  5.  The  higher 
temporal  resolution  allows  analysis  of  more  rapid  phenom- 
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Fig.  7.  Contour  plot  of  the  filtered  image  intensity  at  a  time  of  175  s  that 
identifies  the  distinct  flow  paths  analyzed  in  Figs.  8  and  10. 


Fig.  8.  The  height  of  each  path  defined  in  Fig.  7  vs.  time. 


ena  .  For  example,  the  rapid  rise  of  path  1  is  clearly  resolved 
in  Fig.  8.  Notable  features  of  Fig.  8  include  the  jump  in  the 
height  of  path  2  when  path  1  stops  rising  at  a  time  of  60  s. 
Moreover,  the  connection  between  the  rise  of  path  2  and  the 
fall  of  path  3  at  200  s  is  well  defined.  The  possibility  of  this 
contrast  being  an  artifact  of  the  imaging  can  be  eliminated. 
This  conclusion  is  made  because  of  the  negligible  effect  of 
path  2’s  quick  rise  at  200  s  on  the  other  paths.  The  dominant 
influence  of  path  2  on  path  3  can  be  attributed  to  their  close 
proximity  and  apparent  high  level  of  connectivity. 

For  improved  analysis  of  path  connectivity,  the  gradient 
of  the  fluid  surface  height  with  respect  to  time  is  presented. 
Fig.  9  presents  a  contour  plot  of  the  vertical  velocity  of  the 
fluid  surface  at  194s.  This  point  in  time  corresponds  to  the 
interdependent  rise  and  fall  of  paths  2  and  3,  respectively. 
This  plot  additionally  resolves  any  concerns  that  the  contrast 
of  these  two  paths  might  be  an  imaging  artifact  by  clearly 

x  io~7 


Fig.  9.  Contour  of  the  liquid  surface  vertical  velocity  at  t  =  194  s. 
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Fig.  10.  The  vertical  velocity  of  each  path  defined  in  Fig.  7  vs.  time. 


showing  localized  areas  of  positive  (path  2)  and  negative  (path 
3)  velocity. 

Plots  of  the  fluid  interface  vertical  velocity  in  each  dis¬ 
tinct  path  are  shown  in  Fig.  10.  Though  calculating  gradients 
amplifies  errors  from  signal  noise  and  processing,  the  flow 
features  shown  in  Fig.  8  remain  evident.  A  spike  in  the  veloc¬ 
ity  of  path  2  when  path  1  reaches  the  GDL  surface  at  a  time 
of  60  s  is  present.  Also,  the  sharp  drop  in  velocity  of  path  3 
corresponding  to  the  increased  velocity  of  path  2  is  clear. 


4.  Further  discussion:  transport  mechanism 


From  the  preceding  analysis  of  the  acquired  visualiza¬ 
tion  data,  the  mechanism  for  water  transport  through  the  gas 
diffusion  layers  of  PEM  fuel  cells  is  inferred.  Fig.  1 1  is  an 
illustration  of  the  liquid  interface  generated  when  the  liquid 
surface  passes  through  the  cross-section  generated  by  four 
intersecting  hydrophobic  fibers.  The  capillary  pressure  (Pc), 
which  is  the  jump  in  pressure  across  the  gas/liquid  interface 
due  the  fluid  surface  curvature  (a:),  is  approximated  by: 

Pc  =  ok  (3) 

where 

k  —  Vn  (4) 


where  n  is  the  surface  normal  and  o  is  the  interfacial  tension. 
For  a  square  intersection  of  four  fibers,  the  curvature  can  be 
approximated  with  the  contact  angle  (0),  fiber  spacing  ( S )  and 
fiber  diameter  (D): 


2  a  cos(0) 
S- D 


Assuming  well-dispersed  PTFE  coating  (uniform  0),  the 
transport  within  the  GDL  is  a  process  of  pressure  buildup 
and  breakthrough.  This  process  results  from  the  main  re¬ 
sistance  to  the  progression  of  the  interface  being  capillary 


Fig.  11.  Illustration  of  a  liquid  interface  in  the  cross-section  of  four  inter¬ 
secting  fibers. 


forces  in  the  hydrophobic  structure.  As  water  production  in 
the  catalyst  layer  increases  the  pressure  in  the  liquid,  the  fluid 
surface  slowly  passes  through  the  constrictions  of  the  cross- 
sections  generated  by  intersecting  fibers.  When  the  water  has 
passed  halfway  through  each  constriction,  the  curvature  of 
the  surface  decreases  due  to  the  expanding  radius  of  the  sur¬ 
face.  Thus,  the  liquid  pressure  at  the  interface  drops  and  in¬ 
duces  hydrodynamic  flow  in  the  direction  of  the  expanding 
surface.  Subsequently,  the  liquid  expands  rapidly  above  this 
cross-section  until  the  fluid  interface  contacts  the  next  set 
of  fiber  intersections.  The  process  then  begins  again.  The 
fluid  will  preferentially  pass  through  the  fiber  cross-sections 
featuring  the  greatest  spacing,  as  this  reduces  capillary  pres¬ 
sure  resistance.  An  idealized  schematic  of  the  flow  mech¬ 
anism  is  illustrated  for  the  entire  thickness  of  a  GDL  in 
Fig.  12. 

The  transport  mechanism  depicted  in  Fig.  12  is  consis¬ 
tent  with  the  experimental  observations  presented  earlier  on 
the  evolution  of  water  in  the  various  paths.  The  mechanism 
presented  here  is  different  from  that  proposed  by  Nam  and 
Kaviany  [2]  and  Pasaogullari  and  Wang  [3].  Both  of  these 
works  presented  the  flow  as  an  “upside-down  tree”  capillary 
network.  These  schematics  feature  a  high  number  of  small 
capillaries  dispersed  evenly  within  the  GDL  that  converge 
into  larger  capillaries  and  eventually  result  in  one  very  large 
capillary  that  reaches  the  surface.  However,  the  observations 
herein  promote  the  belief  that  the  transport  is  dominated  by 
fingering  and  channeling  and  features  numerous  “dead  ends”, 
where  water  transport  recedes  when  an  adjacent  breakthrough 
channel  forms.  The  visualization  supports  this  hypothesis  be¬ 
cause  the  dominant  path  is  no  larger  than  an  average-sized 
fiber  cross-section.  In  addition,  after  the  liquid  has  broken 
through  a  constriction  it  is  seen  to  fill  the  entire  void  space 


104 


S.  Litster  et  al.  /  Journal  of  Power  Sources  154  (2006)  95-105 


Fig.  12.  Idealized  schematic  of  the  mechanism  for  liquid  water  transport  in 
PEM  fuel  cell  gas  diffusion  layers. 

surrounded  by  other  fiber  cross-sections.  Thus,  the  conver¬ 
gence  of  small  capillaries  into  successively  larger  ones  is 
not  the  observed  mechanism  in  the  non- woven  gas  diffusion 
layers. 


5.  Conclusion 

Fluorescence  microscopy  has  been  employed  with  a  novel 
methodology  to  provide  new  insight  into  the  dynamic  be¬ 
havior  and  distribution  of  liquid  water  as  it  is  transported 
through  the  gas  diffusion  layers  in  PEM  fuel  cells.  The  new 
quantitative  visualization  methodology  consists  of  pumping 
fluorescein  dye  through  a  GDL  and  correlating  the  light  in¬ 
tensity  emitted  by  the  fluorescein  dye  to  the  height  of  the  fluid 
surface.  The  visualizations  obtained  with  image  processing 
enabled  an  analysis  of  the  connectivity  and  interdependence 
of  the  liquid  motion  in  several  distinct  flow  paths.  The  results 
of  the  analysis  prompted  the  author  to  revisit  the  previously 
adopted  transport  mechanism  hypothesis.  The  experimental 
observations  led  to  the  postulation  of  the  primary  mechanism 
for  liquid  water  transport  in  hydrophobic  GDLs.  The  new  hy¬ 
pothesis  should  prove  valuable  in  the  development  of  models 
for  predicting  liquid  water  transport  in  operating  fuel  cells. 
Moreover,  the  novel  experimental  methodology  provides  an 


additional  method  for  further  characterization  of  gas  diffu¬ 
sion  layers. 
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